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ABSTRACT 

We put limits on the time variation of the electron mass in the early universe 
using observational primordial abundances of D, 4 He and 7 Li, recent data from 
the Cosmic Microwave Background and the 2dFGRS power spectrum. Further- 
more, we use these constraints together with other astronomical and geophysical 
bounds from the late universe to test Barrow-Magueijo's model for the variation 
in m e . From our analysis we obtain —0.615 < Guj/c 4 < —0.045 (3a interval) in 
disagreement with the result obtained in the original paper. 
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Subject headings: primordial nucleosynthesis, cmb, varying fundamental con- 
stants 



Introduction 



Time variation of fundamental constants over cosmological time scales is a prediction 



of th eories that attempt to unify all fundamental interactions 



ries ( IWu and Wang 



1994 



2003 



1983 



Damour et al 



Brax et al 



98 



6: 



2002a. 



Vlaedal Il988l ; iBarr and Mohapatra 



bj), related brane- world t heories ( Youm 



ike s t ring derived field theo- 



1988 



Damour and Polvakov 



200 



2001a 



and Kaluza-Klein theories ( [Kaluza 1921 ; Klein 



b; 



Palma et al. 



1926 



Weinberg 



Gleiser and Taylor 1985 ; Overduin and Wesson 19971 ). In order to study the possi 



ble variation in the fine structure constant or the electron mass, theoretical frameworks 
based on first principles, were developed by different authors ( iBekensteinl Il982j ; iBekenstein 
20021 ; iBarrow et allbooi barrow and Magueiidl2~00~5l ). 



The predicted time behaviour of the fundamental constants depends on which ver- 
sion of the theories is considered. Thus, bounds obtained from astronomical and geophys- 
i cal data are an impo rtant tool to test the validity of these theories. In a previous work 
( IMosquera et al.l 120071 ). we have analyzed the variation in the fine structure constant in the 
context of Bekenstein model. In this paper, instead, we study the variation in the e lectron 
mass (m e ) in the context of the Barrow-Magueijo model (IBarrow and Magueijoll2005l ). Note 
that m e is not a fundamental constant in the same sense as the fine structure constant is. 
Hence, it could be argued that constraints on the time variation of the Higgs vacuum expec- 
tation value (< v >), rather than m e , are more relevant. Moreover, the possibility of a time 
variation of the vacuum expectation value of a field seems more plausible than the time vari- 
ation of a gauge coupling constant. However, in the context of the Barrow-Magueijo model, 
the relevant fundamental constant is m e and thus we will focus on its possible variation. 
The joint variation in the fine structure constant and < v > in the early universe will be 
analyzed in a forthcoming paper. 

Constraints on m e variation over cosmological time scales are available from astronom- 
ical and local methods. The latter ones include geophysical methods (analysis of natural 
long-lived (5 decayers in geological minerals and meteorites) and laboratory measurements 
(comparisons of different transitions in atomic clocks). The astronomical methods are based 
mainly in the analysis of spectra from high redshift quasar absorption systems. Bounds 
on the variation in m e in the early universe can be obtained using data from the Cosmic 
Microwave Background (CMB) radiation and from the abundances of light elements. These 
bounds are not as stringent as the mentioned above but they are important because they 



- 3- 



refer to a different cosmological epoch. 

In this paper, we perform a careful study of the time variation of m e in the early universe. 
First, we use all available abundances of D, 4 He and 7 Li, the latest data from the CMB 
and the 2dFGRS power spectrum to put bounds on the variation in m e without assuming 
any theoretical model. Afterward, we use these bounds and others from astronomical and 
geophysical data, to test Barrow-Magueijo theory. 

In section [2J we use the abundances of the light elements to put bounds on r^ m f , where 
(m e ) is the present value of m e , allowing the baryon to photon density tjb to vary. In section 
[3j we use the three year WMAP data, other CMB experiments and the power spectrum of 
the final 2dFGRS to put bounds on the variation in m e during recombination, allowing also 
other cosmological parameters to vary. In sections HI [5] and [6] we describe the astronomical 
and local data from the late universe. In section we describe the Barrow-Magueijo model, 
and obtain solutions for the scalar field that drives the variation in m e , for the early and 
late universe. In section [8] we show our results. Finally, in section [9] we discuss the results 
and summarize our conclusions. 



2. Bounds from BBN 

Big Bang Nucleosynthesis (BBN) is one of the most important tools to study the early 
universe. The baryon to photon ratio 7]b or equivalently the baryon asymmetry T]b = nj3 ~ nB 
can be determined by comparison between theoretical calculations and observations of the 
abundances of light elements. An independent meth od for determining r i p is provided by data 



from the Cosmic Microwave Background (CMB) (jSpergel et al.l 120031 . 120071 : ISanchez et al. 



20061 ). Considering the baryon density from WMAP results, the predicted abundances are 
highly consistent with the observed D but not with all 4 He and 7 Li. Such discrepancy 
is usually ascribed to non reported systematic errors in the observations of 4 He and 7 Li. 
However, if the systematic errors of 4 He and 7 Li are correctly estimated, we may have 
insight into new physics beyond the minimal BBN model. 

In the currently most popular particle physics models, the lepton and baryon numbers 
are comparable. In this case, any asymmetry between neutrinos and antineutrinos will 
not have a noticeable effect on the predictions of BBN. However, observational data do 
not imply that the lepton asymmetry should be connected to the 'tiny' baryon asymmetry 
T)b- Moreover, a small asymmetry between electron type neutrinos and antineutrinos can 
have a significant impact on BBN since the u e affect the inter-conversion of neutrons to 
protons changing the equilibrium neutron-to-proton ratio from (n/p)° eq = e t~ to (n/p) eq = 
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{ n /p)eq e ~^ e where £ e is the ratio of the neutrino chemical potential to the temperature. 
Consequently, the 4 He abundance changes. In contrast, the D abundance is insensitive to 
£ e 7^ 0. Consist ent with the BBN and CMB data, value s of £ e in the range —0.1 < £ e < 0.3 
are permitted (IBarger et al.ll2003l : ISteigmanll2005l . 120061 ). In this work, however, we assume 
£ e ~ and attribute the discrepancies described above to time- variation of m e or < v >. 

We considered available observat i onal d ata on D, 4 He and 7 L i . For D, we u s ed th e 
values reported bylPettini and Bowenl (120011) ; lO'Meara et al.l (2001); iKirkman et al.l (120 03); 



Buries and Tvtleri Jl998al ]bl): ICrighton et all (120041 ): lO'Meara et ali (hoOfih: 



01 



(2006) . For 4 He, the available observations are re ported by 

( 2007 ) . For 7 Li, we conside r ed the results fromlRyan et al. 



Peimbert et al. 



(2007); 



veira et al. 



Izotov et al. 



(2000); Bonifacio et al. (1997) 



Bonifacio and Molarol (|l997h:lBonifacio et al.l (hood ^ lAsplund et al.l Ki k iBoesgaard et al.l 



(120051 ); iBonifacio et al.l (120071 ) . For the discussion ab out the consistency data check, we refer 
the reader to an earlier work (IMosquera et al.ll2007l ). 



We modified numerical code of Kawano (lKawanolll988l ; lKawanolll992l ) in order to allow 



m e . to vary. The code was also updated with the reaction rates reported by iBergstrom et al. 



( 119991 ) . The main effects of the possible variation in m e in the physics of the first three 
minutes of the universe are changes in the weak rates, in the sum of electron and positron 
energy densities, in the sum of electron and positron pressures, and in the difference of 
electron and positron number densities (see appendix |A] for details). If m e takes a lower 
value than the present one, the primordial abundances are higher than the standards. The 
change is more important for 4 He and 7 Li abundances, where a variation of 10% in m e leads 
to a change of 7.4% and 8.5% in the abundances, while the effect on the D abundance is tiny 
(1.5%). 



We computed the light nuclei abundances for different values of r\ B and 



(m e )o 



and per- 



formed the statistical analysis to obtain the best fit values for these parameters. There is 
no good fit for the whole data set even for , Am { ^ 0. However, reasonable fits can be found 
excluding one group of data at each time (see tabled]). Figures [1] and [2] show the confidence 
contours and 1 dimensional Likelihoods for different groups of data. We obtained that for 
D + 4 He the value of t)b is coincident with WMAP estimation and there is no variation in 
m e within 3a. Moreover, the other groups of data prefer values far from WMAP estimation, 
and for D + 7 Li, the result is consistent with variation in m e within 6a. 

As pointed out in the introduction, in the standard model, D, 4 He and 7 Li abundances 
considered separately predict very different values for the baryon density. Therefore, when 
the three abundances are fitted together, an intermediate value of r\ B is obtained, but the 
value of x 2 is too high. Only when two abundances are considered, we obtain a reasonable 
fit. Furthermore, a high variation in m e which affects mostly the 7 Li abundance is needed 
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Table 1: Best fit parameter values and la errors for the BBN constraints on T^fy- and tjb 
(in units of 1(T 10 ). 





r] B ±a 


, A ™ e ±ct 
(m e ) n 


- 

JV-2 


D + 4 He + 7 Li 


A OQ7+0.047 
-0.097 




9.33 


4 He + 7 Li 


3-6481H 




1.00 


D + 7 Li 


t- oqq+0.287 


u.uuo_ 045 


1.01 


D + 4 He 


U - OOJ -0.355 


-0.022 ±0.009 


1.01 




Fig. 1. — la, 2a and 3a likelihood contours for ; { vs t\b and 1 dimensional Likelihood 
using 4 He + 7 Li data (left) and all data (right) 



to fit D and 7 Li together. On the other hand, D and 4 He are marginally consistent with 
WMAP estimation and therefore no variation in m e is needed to fit both data at the same 
time. Finally, in order to fit the abundances of 4 He and 7 Li, a variation in m e is needed since 
both quantities are affected when m e is allowed to vary. 

As mentioned in the introduction, in this paper we limit ourselves to the context of 
the Barrow-Magueijo model of a varying m e . However, in more general classes of theories 
(Kaluza-Klein, Strings, GUTs, etc), the underlying fundamental constant is the Higgs vac- 
uum expectation value. The dependence of the primordial a bunda nces on the Higgs vacuum 
expectation value has been analyzed by lYoo and Scherrerl t 



have been performed by some of us in earlier works (IChamoun et al.ll2007l ). Besides changes 



200 



^) . Semi- a nalytical analysis 



in m e , a possible variation in < v > modifies the values of the following quantities: the Fermi 
constant Gp, the neutron-proton mass difference Am np , and the deuterium binding energy 
€d ■ The dependence of these quantities with < v > have been described in an earlier work 



n B [10 10 ] AnVK) 



%[io- 10 ] 



Am e /(m e ) 



Fig. 2. — la, 2a and 3a likelihood contours for , Am , e vs 77^ and 1 dimensional Likelihood 
using D + 7 Li data (left) and D + 4 He (right) 



( jChamoun et al.ll2007l ) (see appendix [X] for details). We modified numerical code of Kawano 
in order to allow < v > to vary. The abundances of the primordial elements are much higher 
than the standard value if the Higgs vacuum expectation value during BBN is larger than 
the current value. A variation of 10% in < v > leads to a change of 45%, 25% and 29% in 
the 4 He, 7 Li and D abundances respectively. Since D is a residual of 4 He production, a great 
change in 4 He also leads to an important change in D. The changes in the abundances are 
greater than in the case where only m e is allowed to vary. 

In the case of < v > , we performed the same analysis for the same groups of data we 
considered for m e . As in the case of m e , variation there is no good fit for the whole set of 
data. However, reasonable fits can be found excluding one group of data at each time (see 
table [2]). Figure [3] shows the confidence contours and 1 dimensional Likelihoods for different 
groups of data. 



Table 2: Best fit parameter values and la errors for the BBN constraints on and 
(in units of 10~ 10 ). 





Vb ± CF 


A < v > ± a 

<v>n 


■ 

N-2 


D + 4 He + 7 Li 


4.275 ±0.097 


0.006 ±0.002 


9.27 


4 He + 7 Li 


o *roq+0.132 
°- ' ZO -0.124 


0.008 ±0.001 


1.00 


D + 7 Li 


K 1 QO+0.242 
0.lOJ_ Q 231 


n i qc+0.015 
u.±oo_ 009 


1.01 


D + 4 He 


6.324^| 74 


0.004 ±0.002 


1.04 
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We obtain that for D + 4 He the value of tjb is consistent with WMAP estimation and 
there is no variation in < v > within 3a. Moreover, the other groups of data prefer values not 
consistent with WMAP results. For D + 7 Li, the result is consistent with variation in < v > 
within 6a. The results are similar to those obtained in the case where m e is the varying 
constant: i) no reasonable fit for the three abundances; ii) D and 4 He can be well fitted with 
null < v > variation; iii) D and 7 Li need a huge variation in order to obtain a reasonable fit. 
However, the bounds on variation in < v > are more stringent than the bounds obtained 
when only m e was allowed to vary (see tabled]). This could be explained since variations in 
< v > lead to greater changes in the theoretical abundances than variation in m e . 




t 1b [10" 10 ] A<v>/<v> 1Tb [10" 1d ] A<v>/<v> 

Fig. 3. — lcr, 2a and 3a likelihood contours for ^ < ^ l> vs ttr and 1 dimensional Likelihood 
using D + 7 Li (left) and D + 4 He (right) 



3. Bounds from CMB 

Cosmic Microwave Background (CMB) radiation provides valuable information about 
the physical conditions of the universe just before decoupling of matter and radiation, and 
thanks to its dependence upon cosmological parameters, it allows their estimation. Any 
change in the value of m e affects the physics during recombination, mainly the redshift 
of this epoch, due to a shift in the energy levels and in particular, the binding energy of 
hydrogen. The Thompson scattering cross section, which is proportional to m~ 2 , is also 
changed for all particles. Therefore, the CMB power spectrum is modified by a change in 
the relative amplitudes of the Doppler peaks, and shifts in their positions (see appendix 
IB1 for details). Changes in the cosmological parameters produce similar effects. In the 
recombination scenario, the only effect of varying < v > is a change in the value of m e . 
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Pr evious analysis of the CMB data including; a possible variation in m e have been performed 

heWMAP 
2007h . and 



by lYoo and Scherrerl (120031 ) ; iKujat and Scherrerl (120001 ) . In this p aper, we use t 



3-year temperature and temperature- polarization power sp ectrum (ISpergel et al 



other CMB exper iments such as CBI (IReadhead et a 



BOOMERANG (Piacentini et al.ll2006l: Ijones et al. 



20041 ). ACBAR (IKuo et al 



2004J), and 



20061 ). and the power spectrum of the 



2dFGRS (jCole et al.l 120051 ) . We consider a spatially-flat cosmological model with adiabatic 
density fluctuations. The parameters of our model are: 



P = {Q B h 2 , QcDAih 2 , 6, r r , 



Am P 



m 



where flcDAih 2 is the dark matter density in units of the critical density, O gives the ratio 
of the comoving sound horizon at decoupling to the angular diameter distance to the surface 
of last scattering, r re is the reionization optical depth, n s the scalar spectral index and A s is 
the amplitude of the density fluctuations. 

We use a Markov Chain Monte Carlo method to explore the parameter space because the 
exploration of a multidimensional parameter space with a g rid of points is computa tionally 
prohibitive. We use the public available Cosmo MC code of iLewis and Bridle! ( 120021 ) which 
uses CAMB (Lewis et al.l boooh and RECFAST JSeaeer et al.l Il999h to compute the CMB 
power spectra, and we have modified them in order to include the possible variation in 
m e at recombinati o n. W e ran eight different chains. We used the convergence criterion of 



Raftery and Lewis! (119921 ) to stop the chains when R — 1 < 0.0044 (more stringent than the 
usually adopted value). Results are shown in table Eland figure HI Figure 0] shows a strong 
degeneracy between m e and G, which is directly related to H , and also between m e and 
Qsh 2 , and m e and flcDAdh 2 . 

We have also performed the analysis considering only CMB data. The strong degeneracy 
between m e and H made the chains cover all the wide H prior, making it impossible to 



Table 3: Mean values and errors for the main and derived parameters including m e variation 
(Hq is in units of km s _1 Mpc -1 ). 



Parameter 


Mean value and la error 


n B h 2 


0.0217 ±0.0010 


e 


1.020 ±0.025 


Am e 
(m e ) n 


-0.029 ±0.034 


A s 


3.020 ±0.064 



Parameter 


Mean value and la error 


&<CDAlh 2 


1 00fi +0 0085 
u.iuuu_ 0086 




n nqi +0.013 
u.uyi„ 014 


n s 


0.960 ±0.015 


H 
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Fig. 4. — Marginalized posterior distributions obtained with CMB data, including the 
WMAP 3-year data release plus 2dFGRS power spectrum. The diagonal shows the pos- 
terior distributions for individual parameters, the other panels shows the 2D contours for 
pairs of parameters, marginalizing over the others. 



find reliable mean values and errors. Hence, we added a gaussian prior to Hn , which was 
obtained from the Hubble Space Telescope Key Project (IFreedman et al.l 120011 ). and chose 
the values of the mean and errors as those inferred from the closest objects in that paper, 
so we could neglect any possible difference between the value of m e at that redshift and the 
present value. In this way, we post-processed the chains and found limits that are consistent 
with those of the first analysis, revealing the robustness of these bounds. However, the most 
stringent constraints were obtained in the first analysis (see figure [5]). 



Finally, we comment that Planck will be the fir st mission to map the entire CM B sky 
with mJy sensitivity and resolution better than 10' (IThe Planck Collaboration! 120061 ). Such 
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Fig. 5. — la - , and 2a confidence levels contours obtained with CMB data with and without 
data of the 2dFGRS power spectrum. 

resolution will allow to see into the damping tail of the anisotropy spectrum, arou nd the 



third and fourth peaks, with a precision almost only limited by cosmic variance (IWhite 



20061 ). This will enable a very precise estimation of the baryon-to-photon ratio from the 



relative height of the peaks in the spectrum. 



4. Bounds from Quasar Absorption Systems 



Quasar absorption systems present ideal laboratories to search for any temporal varia- 
tion of the fundamental constants over cosmological ti me scales. In particular, a method fo r 
constraining the variation in /i = ^ was developed by IVarshalovich and Levshakovl (119931 ). 
It is based on the fact that wavelengths of electron-vibro-rotational lines depend on the re- 
duced mass of the molecules, with different dependence for different transitions. In such way, 
it is possible to distinguish the cosmological redshift of a line from the shift caused by a vari- 
ation in fi. The rest-frame laboratory wavelength, A°, can be related to those in the quasar 
absorption system, A i; as |fr = (1 + z a b s )(l +K i —^), where z a b s is the absorption redshift and 
Ki is the coefficient which determines the sensitivity of the wavelength Aj. Using observations 
from H 2 absorption syste ms at high redshift and laboratory measurements, several authors 



obtai ned constraints on fi fjPotekhin et al.lll998t iLevshakov et al.ll2002l ; llvanchik et al.ll2003 



20051 ). The most up-to-date available measurements for each redshift are listed in table 0] 
and will be considered to test Barrow-Magueijo model. 



Another method for constraining variation in fundamental constants is based on the 
comparison between the hyperfine 21 cm absorption transition of neutral hydrogen (y a ) and 
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Table 4: The table shows the absorption redshift, the value of with its corresponding 
error (in units of 10~ 5 ), and the reference obtained comparing molecular and laboratory 
wavelengths. 



Redshift 


, A ™ e ±a 

(m e )n 


Reference 


2.8 


-6.25 ± 13.70 


Potekhin et al. f 19981 


3.02 


-1.40 ±0.83 


Ivanchik et al. (2005) 


2.6 


-2.11 ± 1.39 


Ivanchik et al. (2005) 



an optical resonance transition (uh). The ra tio ^ is proportional to x — a 2 g p ^ where g p is 
the proton g factor (ITzanavaris et al.ll20071 ). Thus, a change of this quantity will result in a 



difference in the redshift measured from 21 c m and optical absorpti o n line s as — — Zopt 



(l+22l) ' 

Since we are working in the context of the iBarrow and Magueijol (j2005l ) model, the only 
fu ndamental constant whi ch is allowed to vary is m e . Table [5] shows the bounds obtained 
by ITzanavaris et al.l (120071 ) combining the measurements of optical and radio redshift. This 
method has the inconvenience that it is difficult to determine if both radio and optical lines 
were originated at the same absorption system. Thus, a difference in the velocity of the 
absorption clouds could hide a variation in x. 



Table 5: The table shows the absorption redshift, the value of j^y with its corresponding 
error (in units of 10~ 5 ), obtained comparing radio and molecular redshifts (ITzanavaris et al. 



2007J). 



Redshift 


(m e ) n 


0.24 


1.21 ± 2.10 


0.31 


-0.61 ±4.27 


0.40 


3.22 ±3.15 


0.52 


-2.95 ± 1.05 


0.52 


0.26 ±3.67 



Redshift 


(m e ) n 


1.78 


-2.59 ±0.90 


1.94 


3.30 ± 0.44 


2.04 


5.20 ±2.76 


2.35 


-2.54 ± 1.82 



5. Bounds from Geophysical Data 

The half-life of long-lived j3 decayers has been determined either in laboratory measure- 
ments or by comparison with the age of meteorites, as found from a decay radioactivity 
analysis. The most stringent bound on the variation in the half life, A, proceeds from the 
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comparison of 187 Re dec ay in the Solar System forma tion and the present (lOlive et al.ll2004j ): 
^ = (— 0.016 ± 0.016). ISisterna and Vucetichl (1990) derived a relation between the shift in 
the half-life of long lived (3 decayers and a possible variation between the values of the fun- 
damental constants values now. In this paper, we only consider m e variation and therefore 



AA 
A 



, Am e 

'(m e ) 



where a = —600 for 187 Re. 



6. Bounds from Laboratory 



Comparison between frequencies of different atomic transitions over time are useful tools 
to put stringent bounds on the variation in fundamental constants at present. In particular, 
the hyperfme frequency of cesium can be approximated by z/ Cs ~ g Cs ^a 2 R y F Cs (a) (where 
gcs is the nuclear g factor, R y is the Rydberg constant expressed as a frequency and F Cs (a) 
is a dimensionless function of a and does not depend on m e at least at first or der), while 
optical transition frequen c ies can be expres sed as v opt — R y F(a). Several authors (iBize et al. 



20031 ; iFischer et al.l 12004 ; iPeik et al.l 12004 ) have measured different optical transitions and 
compared them with the frequency of the ground state hyperfme splitting in neutral 133 Cs. 



These measurements can be used to constrain the variation in 
different experiments are listed in table El 



(m e ) • 



Constraints from 



Table 6: The table shows the compared clocks, the value of j^y with its corresponding error 
fin units of 10~ 15 ), the time interval for which the variation was measured and the reference. 



Frequencies 


,™\ ± a [10- 15 yr _1 l 


At[yr] 


Reference 


Hg + and Cs 


0.2 ±7.0 


5 


Fischer et al. (2004) 


Yb+ and Cs 


1.2 ±4.4 


2.8 


Peik et al. (2004) 


Hg + and Cs 


0±7 


2 


Bize et al. (2003) 



The Model 



We now analyze the Barrow-Magueijo model for the variation in m e . We solve the 
equation of the scalar field (0) that drives the variation in m e in this model. We consider 
that the variations in <ft are small and they do not pr oduce significants cont ributions to 
the Friedmann equation. As we did in a previous work (IMosquera et al.l 120071 ) . we build a 
piecewise approximate solution by joining solutions obtained by keeping only some terms 
of the Friedmann equation, relevant in the following domination regimes: a) radiation and 
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matter, and b) matter and cosmological constant. In such way, solution a) can be applied to 
nucleosynthesis and recombination of primordial hydrogen whereas solution b) is appropriate 
for quasar absorption systems, geophysical data and atomic clocks. 

Defining the variable •& as di? = dr/a, where r = H t, and t is the cosmic time, the 
expression for the scale factor in the radiation and matter regime is: 

a RM {$) = ^ 2 +v^> (2) 
and the relationship between r and t? is: 

r(tf) = 5^ 3 + ^ (3) 

The solution for the scale factor in the matter and cosmological constant regime can be 
written as: 

2/3 

(4) 



'n xl/3 



3 

sinh ( - V"A (t — r ) + sinh" 1 




where To = H to, and to is the age of the universe. To obtain the last solution we have 
considered that the scale factor must be a continuous and smooth function of time. 

In the Barrow-Magueijo model, m e is controlled by a dilaton field defined by m e = 
(m e ) exp(0), and variations in m e occur no sooner the universe cools down below m e thresh- 
old. The minimal dynamics for <p is set by the kinetic Lagrangian 

= (5) 

where a; is a coupling constant. From this Lagrangian, the equation of motion of the scalar 
field can be derived as: 

(0a 3 ) = -Mexp[<j>] (6) 

with M ~ peQa^c^/uj. This is a second order equation for 0, with the boundary condition 
0o — 0. If the mass variations are small, ~ 1 can be set to obtain an analytical expression 
for 0. 

For convenience, we define (3 = \VL m VLr 1 ^ £ = VtmVtr^^Vf^^ , 7 = Vt^Vb^^ 2 and 

m = MI 

C = sinh-'r'^-ZK) (T) 

Provided = ln(m e / (m e ) ) — Am e / (m e ) , the expressions for the variation in m e in the two 
regimes are: 
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for r < T\ (where T\ is defined by Or^ti) = a Mc( T i) — 
(ExS)): 



1/4 



see 



Mosquera et al. 



Am e 



2 M 1 

3 fi m 



A Q m 



-2 In 



2(/ft9 + i; 



+ 



2 +?/(*)>/i+! 



y/T+ 



T 



7 



1 



1 1 / y/r+e-r 

2 n V /ft? y 2 VvT+I+i; 4/3tf 4(/M+l) 



(e-l)v / r+e + 



3 7 



The relationship between r and is given by Eq.Q. 
for t > T\ 



Am e 
(m e ) 



(r) 
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M 2 
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QaT coth ( C + — y tt&T J — — In 
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C + In f 7 + ^/TT 



7' 



+ 



A 2v^Xv 
i? 3fi m 



coth(<7+~V^ + ^lll! 



(9) 



where A is an integration constant. 



8. Results 

The model described in section [7] predicts the variation in m e as a function of time, 
and has two independent dime nsionless parameters M/Hn and A/Hq. We do not fix A/Hq 



to zero as previous works did (IBarrow and Magueijoll2005l ). To constrain these parameters, 
we use the data described in the previous sections. We perform a \ 2 test to obtain the 
best fit parameters of the model. In order to obtain the parameters consistently with our 
assumption that the energy density of the field <\> (e^ = ^f</> 2 j can be neglected in the 
Friedmann equation, we add to the x 2 expression, a term that controls that the contribution 
of (f) to the Friedmann equation will be less important that the radiation term, right after 
m e threshold is crossed. The result of the statistical analysis shows that there is no good fit 
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for the whole data set. We repeat the analysis excluding one group of data at each time. 
We found that reasonable fits can be obtained excluding the quasar at z — 1.94 of table [5] 
and that the data from nucleosynthesis is crucial to determine the value of A/H . Besides, 
the group of data from table H] is important to determine the value of M/Hq. The results 
are shown in table [71 

Table 7: The table shows the best fit parameters of the model (excluding entry 7 of table [5]). 
The value for the Mj Hq parameter is in units of 10 -6 , and the value for the A/Hq parameter 
is in units of 10~ 13 . 



M/Hl 


A/H 


Glu/c 4 


X 1 - 
7V-2 


-7 sn+ 210 

' - OU -2.02 


3 fiD +1 - 44 
°- DU -1.50 


-n S36+ - 097 

U.OOU_ 093 


1.14 



Since the Barrow-Magueijo model is written in terms of the coupling constant u, we 
derive its best value from the previous constraints. Since M ~ p^a^c^juj we obtain the 
following relationship: 

Gu_3 ( f Uc \ m e (My 1 

where /fie is the fraction of the total number of baryons in the form of He, and can be 
written as a function of the total observed mass abundance of He (Mne/Mn). According to 
the values of M/Hq from table [7J and using / He = 0.19 (taking M He /M H = 0.24), we obtain 
the bounds on the dimensionless quantity ^ presented in table [71 



9. Summary and Conclusions 

In this paper, we have put limits on the time variation in the electron mass at primordial 
nucleosynthesis time using observational primordial abundances of D, 4 He and 7 Li, and we 
have analyzed in detail the consequences of considering different groups of data. We have 
also considered the variations in < v > during BBN and analyzed the differences with the 
variations in m e during the same epoch. Additionally, we have used the three year data from 
the Cosmic Microwave Background and the final 2dFGRS power spectrum to obtain bounds 
on the variation in m e at recombination, and an estimation of the cosmological parameters. 
Together with other bounds on the variation in the late universe, that come from quasar 
absorption systems, half-life of long-lived (3 decayers, and atomic clocks, we put constraints 
on the Barrow-Magueijo model for the variation in m e . We have improved the solutions by 



-16- 



taking into account the detailed evolution of the scale factor and the complete solution for 
the scalar field that drives the variation in m P . 



In the original paper (IBarrow and Magueijoll2005l ) some approximations in the evolution 
of the scale factor are assumed with the consequent simplification in the solution for the 
scalar field. Another improvement of our derivation is that we have not neglected the first 
integration constant, which is the most contributing part in the early universe. In fact, 
integrating Eq. ([()]) once, we can write: 



<> a 3 = -M \ t - — I 



-) 

Ml 



-M (t + 8.47 x 10 2 yr) 



(11) 



where we have used the best fit values for the parameters M/Hq and A/ Ho, and h = 0.73. 
Note that the second term in the right hand side of Eq. fllll) is dominant in the early universe, 
in particular, during nucleosynthesis. 



Barrow and Magueijol (120051 ) presented a bound of G\u\ > 0.2 (with c = 1). They 
obtained such constraint using bounds from quasars at z ~ 1, whereas we use all the available 
bounds on the variation in m e at different cosmological times. In appendix [C] we briefly 
discuss the difference in both analysis. From data supporting the weak equivalence principle, 
they obtain G\u\ > 10 3 while we obtain —0.615 < Guj/c 4 < —0.045 (3a interval) using data 
from different cosmological time scales. More research both on time variation data and on 
the bound from WEP is needed to understand this discrepancy. 

Finally, we remark that, at 2cr, the value of us is negative. This should not be surprising. 
Indeed, negative kinetic terms in the Lagra ngian have alre ady been considered in k-essence 
models with a phantom energy component (ICaldwellll2002l ). 
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A. Physics at BBN 



In this appendix we discuss the dependences of the physical quantities involved in the 
calculation of the abundances of the light elements with m e and variation in the Higgs vacuum 
expectation value. We also discuss how this quantities are changed within the Kawano Code. 



A.l. Variation in the electron mass 



A change in the value of m e at the time of primordial nucleosynthesis with respect to its 
present value affects derived physical quantities such as the sum of the electron and positron 
energy densities, the sum of the electron and positron pressures and the difference of the 
electron and positron number densities. In Kawano's code, these quantities are calculated 
as follows: 



7T 

T 

m e c2 



he 



3 13 



m P cr 



Pe~ + Pe+ 

Pe- + Pe+ 
C 2 

z 3 (n e - - n e +) 



\m e c 



2n4 
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(he) 



^(-l) n+1 ch(n0 e )M(nz) 



\m e c 



2\4 



(he) 



-1 



\n+l 



nz 



-ch (n0 e ) N(nz) 



^(-l) n+1 sh(n0 e )L(nz) 



(Al) 
(A2) 
(A3) 



where z = ^ff^, <p e is the electron chemical potential and L(z), M(z), and N(z) are related 
to the modified Bessel function (jKawanol Il988l ; iKawanol Il992l ). In order to include the 
variation in m e we replace, in all the equations, m e by (m e ) ^1 + ) ■ The change in 



(m e ) 

these quantities, due to a change in m e , affects their derivatives and the expansion rate 
through the Friedmann equation. The n <-» p reaction rates and the other weak decay rates 
are changed if m e varies with time. The total n — > p reaction rate is calculated by: 



A = K 



dE e E ePe (E e + Am. 



np ) 



+ K 



dE e E e p e (E e - Am 



np! 



(1 _|_ e Ee/Ty} (I _|_ e -(E e +Am np )/T v -^ ^ ^ _|_ g-S e /T 7 ) + e (E e -Am np )/T u -£ l ^ 

where E e and p e are the electron energy and momentum respectively, Am np is the neutron- 
proton mass difererences, K is a normalization constant proportional to G 2 Fl and E e = 
(p 2 e + m 2 e ) 1 12 . 

It is worth to mention that the most important changes in the primordial abundances 
(due to a change in m e ) arrives from the ch ange in the weak rates rather than the change in 
the expansion rate (jYoo and Scherrerl 120031 ) . 



(A4) 
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A.2. Variation in the Higgs vacuum expectation value 



If the the value of < v > during BBN is different than its present value, the electron mass, 
the Fermi constant, the neutron-proton mass difference and the deuterium binding energy 
take different values than the current ones. The electron mass is proportio nal to the Higgs 
vacuu m expectation value. The Fermi constant is proportional to < v >~ 2 fjDixit and Sher 
19881 ). This dependence affe ct the n «-» p reaction ra tes since K ~ G 2 F . The neutron-proton 
mass difference changes by (IChristiansen et al.lll99ll ) 



5 Am 



rip 



Am. 



1.587- 



up 



.A < v > 

<v> 



(A5) 



affecting n <-> p reaction rates (see Eq. (lA4l) ) and the initial neutrons and protons abundances: 

1 1 



Y„ 



\ _|_ gAtn„p/T 9 +$ 



I _|_ e -Am„ p /r 9 -^ 



(A6) 



where Tg is the temperature in units of 10 9 K. In order to include these effects we replace 

The deuterium binding energy must also be corrected by 



Am np by Am Ti 

Aeo _ a<v> w ] iere j s a m odel dependent constant. In this work we assume, following 

— £t> <V>0 1 1 1 r ' ° 

Chamoun et al.l (120071 ). k = —0.045. This correction affects the initial value of the deuterium 



abundance. Once again we replace by €d fl + \ in the code. 



B. Physics at recombination epoch 

During recombination epoch, the ionization fraction, x e , is determined by the balance 
between photoionization and recombination. The recombination equation is 

where 

C = (! + KK 2sM n ls ) 

(l + K((3 c + A 2sM )n ls ) { ) 

is the Peebles factor, which inhibits the recombination rate due to the presence of Lyman- 
a photons. The redshift of these photons is K = with X a = and A 2s ,i s is 

the rate of decay of the 2s excited state to the ground state via 2-photon emission, and 
scales as m e . Recombination directly to the ground state is strongly inhibited, so the case B 
recombination takes place. The case B recombination coefficient a c is proportional to m e 3 ^ 2 . 
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The photoionization coefficient depends on a c , but it also has an additional dependence on 

^l-nnipkTr, 

The most important effects of a change in m e during recombination would be due to its 

Thomson 
hydrogen B\ = 1 a 2 m e c 2 . 



3/2 

^ -■^ , e -B 2/fc T m (B3) 

n e during 

influence upon Thomson scattering cross section <jt = 8 Z h2 f a 2 , and the binding energy of 

o in 



C. Different limits on Guo 



In this appendix we compare the limits obtained by lBarrow and Magueijol (120051 ) on Guo 
with our bounds. We stress that we have performed a \ 2 using all available observational 
and experimental data, while iBarrow and Magueijol ( 20051 ) consider | — | < 10~ 5 for data at 
redshift of order 1. Moreover, most exact individual bounds from quasar absorption systems 
are not consistent with null variation at least at la. 

Let us consider for example t he last entry of table [1 1 —3.5 x 10~ 5 < =^ < —0.72 x 10~ 5 . 
Using the same approximation as Barrow and Magueijol ( 2005 ). we find that —0.28 < Glu < 
—0.05 for this measurement which is of the same order of magnitude as obtained considering 
all data. 
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